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A Study  of  the  Error  of  Discretization  in  the  Air  Force 
Global  Weather  Central  Boundary  Layer  Model 


I.  INTRODUCTION 

This  report  summarizes  the  results  of  an  investigation  that  was  undertaken  to 
examine  how  the  method  of  discretization  and  the  structure  of  the  computational 
grid  in  the  vertical  direction  employed  in  the  Air  Force  Global  Weather  Central 
Boundary  Layer  Model  (AFGWC-BLM)  might  influence  the  forecast  accuracy  of  the 
model. 

AFGWC-BLM  simulates  the  evolution  of  the  planetary  boundary  layer  over  a 
continental  area  through  the  numerical  solution  of  an  initial-boundary  value  prob- 
lem posed  on  a set  of  grid  points  distributed  within  the  domain  of  interest.  A 
measure  of  forecast  accuracy  may,  therefore,  be  obtained  from  the  aggregate  of  the 
differences  between  the  predicted  and  observed  values  of  the  variable  of  interest. 
The  effects  that  a particular  part  or  substructure  of  the  model  have  upon  forecast 
performance  can  be  studied  by  regarding  the  whole  model  as  a system  and  bv  car- 
rying out  numerical  experiments  in  which  the  part  under  investigation  is  kept 
under  control.  Such  a system  approach  for  the  study  of  prediction  models  was 
proposed  as  a standard  procedure  to  be  used  for  diagnosis  and  improvement  of 
dynamical  prediction  models  in  a previous  report.  In  the  present  study,  a 
different  approach  is  used  for  the  following  reasons:  AFGWC-BLM  is 
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constructed  with  the  aid  of  many  modeling  assumptions  and  many  mathematical 
approximations.  This  implies  that  the  model  is  endowed  with  as  many  independent 
sources  of  error  as  the  number  of  such  artifices.  To  make  matters  worse,  there 
is  virtually  no  way  of  ascertaining  how  well  these  artifices  represent  the  real  state 
of  the  atmosphere  at  any  given  specific  time  and  place.  It  is,  therefore,  extreme- 
ly difficult  to  accurately  trace  the  forecast  error  of  the  model  to  a particular 
source  or  cause  without  extensive  experimentation. 

However,  it  is  possible  to  make  a reasonable  assessment  of  the  influence  or 
significance  that  a particular  substructure  of  the  model  has  upon  forecast  accuracy, 
if  one  does  not  insist  upon  an  answer  that  establishes  a cause-and-effect  relation- 
ship. This  is  done  by  isolating  the  substructure  of  interest  from  the  model  and  by 
examining  its  characteristics  independent  of  the  rest  of  the  model.  This  method  of 
investigation  may  be  applied  if  evaluation  does  not  require  observational  data  of  the 
real  atmosphere,  and  is,  therefore,  particularly  suitable  for  studying  the  nature 
of  the  mathematical  structure  of  the  model.  The  subject  of  this  investigation  is 
one  example  of  such  a method.  However,  as  will  be  shown  in  Section  2,  even  this 
method  requires  considerable  simplification  and  qualification  in  order  to  obtain 
definite  conclusions. 


2.  SCOPE  AND  METHOD  OF  INVESTIGATION 

To  explain  the  rationale  underlying  the  investigation,  let  us  start  by  consider- 
ing a typical  prediction  problem:  Given  a suitable  set  of  initial  and  boundary  condi- 
tions,  we  search  for  the  solution  y(x,  t)  in  a bounded  continuous  domain  D,  xoR, 

0 't  = T,that  satisfies  differential  equation 

3y  ,/■*  . x c R 

= f(x*  y.  o<t?T  ’ (1> 

in  which  f is  generally  a nonlinear  function  of  the  variables  concerned.  It  is 

assumed  that  the  existence  and  uniqueness  of  the  solution  is  beyond  doubt  and, 

therefore,  presently  out  of  question.  We  shall  denote  by  y]}  the  solution  satisfying 

Eq.  (1)  and  the  other  accompanying  conditions. 

The  present  state  of  the  art  in  dealing  with  this  problem,  except  for  a very  few 

cases  of  extreme  simplicity  and  limited  practical  value,  usually  requires  two  steps: 

discretization  and  linearization.  The  continuous  domain  D is  replaced  by  a finite 
- * — ) 

set  A of  grid  points  (x  , t ) formed  by  nets  covering  D.  The  value  v(x.,  t ) at  each 
’ ^ in 

of  these  grid  points  is  sought.  Function  f is  linearized  so  that  the  discrete  analog 

of  the  differential  equation  is  reduced  into  a system  of  algebraic  equations  of  the 
form 
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(2) 


v.(x. , t , ) = y.(x.,  t ) + F(v.(x,  , t ),  x,  , t ) . 
A i n+ 1 Ain  A k n k n 


The  set  1 v^lx.,  t|i+  m>,  x.:R,  lSmSNj.  constitutes  the  prediction,  and  will  be 
referred  to  as  the  discrete  solution. 

A direct  measure  of  the  accuracy  of  such  a prediction  can  be  readily  obtained 

from  the  differences  v^lx.,  t r ) - yD(x.,  tn+m',  except  for  the  fact  that  there  is 

generallv  no  analvtic  solution  whose  values  ran  be  obtained  as  readily  as  those  of 

v ,(\-  , t ) for  anv  practical  function  f(v,  x,  t).  A substitute  is  therefore  proposed 
A i n+  n i 

that,  while  not  directlv  representative  of  prediction  accuracy,  expresses  the  degree 
of  inaccuracv  incurred  through  the  process  of  discretization  and  can  be  used  for 
assessing  the  merit  of  anv  particular  procedure  in  comparison  with  other  proce- 
dures performing  the  same  function.  The  substitute  is  the  difference  between  two 
steadv-state  solutions  (\\,  and  y-,  respectively)  of  llqs.  (1)  and  (2)  on  the  set  of 
grid  points,  where  \y  and  \v  are  obtained  by  removing  the  time  dependence  from 
both  equations.  The  quantity  measures  one  aspect  of  inaccuracv,  even  for  the  time- 
dependent  problem,  since  it  induces  a change  with  time  in  the  discrete  svstem,  even 
when  the  continuous  system  is  in  the  stationary  state.  We  shall  call  it  the  discreti- 
7.  ition  error  in  the  steadv-state  solution. 

\l  GWC-BI.AI  emplovs  a computational  grid  that  has  a net  of  uniform  spacing 
on  a stereographic  projection  in  tiie  horizontal  direction  and  a constant  time-step 
for  the  entire  period  of  time  integration.  In  the  vertical,  however,  the  grid  is  made 
up  of  eight  unevenly  distributed  levels.  They  are  the  surface  level,  to  which  the 
height  of  1.2?>  m above  ground  level  (AGI.)  is  assigned  in  the  computation,  50  m, 

1 50  m,  iOO  m,  tiOO  m.  !>00  m,  1200  m,  and  1600  m AGI.. 

1 he  model  simulates  the  vertical  turbulent  transfer  of  heat  and  moisture  with 
the  use  of  eddv  liffusivitv.  The  eddy  diffusivitv  is  defined  to  be  dependent  on  the 
eloeitv  nd  temperature  distributions  in  the  vertical  and,  therefore,  varies  with 
hi  tune  i nd  height. 

nt  i ,11  spatial  derivatives  in  the  continuum  are  replaced  by  finite  differ- 
ent . on  the  computational  grid,  our  attention  was  drawn  particularly  to  the  manner 
in  which  discretization  is  carried  out  in  the  vertical.  A simple  model  that  incor- 
porates ill  aspects  of  concern  is  obtained  in  the  Dirichlet  problem  of  heat  conduc- 
tion defined  below  to  find  viz)  that  satisfies  (A)  and  (B): 


(A)  Differential  liquation  \ 

(K<z)v,',  S(z)  I 

1 > z > 0 <31 

K(z)  - Kn  > 0 

with  v , = dy/dz  l 

(B)  Boundary  Conditions  y(0)  = y and  y(l)=  y..  ' 


— 


Here,  two  parameter  functions  K(z)  and  S(z)  are  so  prescribed  that  the  existence  of 
a continu 
given  by 


a continuous  solution  is  insured.  It  is  readily  shown  that  solution  y<.(z)  may  be 


z C i + J S(v)dv 

>s(z)  = yo + i0  — m du  • 


where 


C1  = 


(.Yi  - y0) 


1 / z 
~ ‘^0  ' Jo 


S(u)du/K(z) 


dz 


(4) 


dz/K(z) 


The  solution  can  be  easily  evaluated  when  K(z)  and  S(z)  are  simple  functions  of  z. 

Next,  the  discrete  version  of  the  Dirichlet  problem  (3)  is  obtained  by  using  a 
given  method  of  discretization.  The  discrete  version  takes  the  form  of  a system 
of  equations  whose  solution  isvy(z.),  which  together  with  Vg(z.)  defines  the  discreti- 
zation error  in  the  steady-state  solution.  Comparisons  of  discretization  errors 
among  a class  of  comparable  methods  of  discretization  will  then  enable  us  to 
assess  the  relative  merits  of  individual  methods. 

Because  of  the  limited  significance  of  the  discretization  error  in  the  steady- 
state  solution,  we  did  not  believe  it  worthwhile  to  generalize  the  parameter  func- 
tions or  the  boundary  conditions.  If  the  results  on  simple  parameter  functions 
afford  an  unequivocal  evaluation,  we  feel  that  they  should  be  sufficient  for  us  to 
make  choices.  The  investigation  is  thus  confined  to  cases  that  arise  from  com- 
bining one  eddy  diffusivity  and  one  source  function  from  the  sets  listed  below.  Pro- 
files of  some  of  the  eddy  diffusivities  and  the  source  functions  are  shown  in  F igures 
1 and  2 respectively.  Those  that  are  not  shown  maybe  inferred  easily  from  those 
shown. 

(A)  Eddy  Diffusivity,  K(z) 


Kj  = -49(z  - 
K2  = -49(z  - 
K.j  = -49(z  - 

K.  = 100(z  - 

4 

K.  = 100(z  - 

5 

K„  = 1 00(z  - 
b 


25)  + 28.  5625 

.2 


375) 


20.  140625 


, 50)*  + 13.25 

, 32)2  + 1.00 

,2 


32) 

32P 


0.  10 
0.  01 


(B)  Source  Function,  S(z) 


S1  = 0 

S2  = 10(1  - z2) 

S = -40(z  - 1/2)2  + 


10 


S,  = 20  - S„ 


S5  = 


10z(2 
= -Sn 


z) 


00  0'  02  03  04  05  06  07  08  09  10  00  01  2 03  4 d1-,  6 7 R '9 

YS—  YS  — 

Figure  3.  Sample  Continuous  Functions.  ys>  Identifier!  hv  Pair  of  the 
Parameter  Functions  (K.,  S.) 

i .1 


Two  different  finite-difference  formulations  of  the  left-hand  side  of  the  equation 
of  heat  conduction,  Kq.  (3),  and  three  different  computational  grids  were  considered. 
Comparisons  were  made  to  put  a separate  focus  on  e:  ch  of  two  different  aspects  of 
discretization;  namely,  the  finite-difference  formulation  and  the  structure  of  com- 
putational grids  represented  bv  lensitv  and  distribution.  Faeh  subieet  will  be 
treated  separately  in  Section  3. 


0'l/'  0i7 
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1.  Kl  'I  l.l  S Ol  IW  KSI  l(,  VI  ION 


3.1  linin'  Dilim'iict'  for  imil.il  ioii- 


VVe  consider  the  following  two  finite-difference  formulations  for  the  expres- 
sions (K(zly,>/  on  a computational  grid  (z.  , v 1,2,...,  X).  Tvpe  1 is  suggested  by 
writing  <K(z)y  ) * K(z)(y  1 + kV*.  Thus,  lOq,  (3)  becomes  discretized  as 


2 K. 

i 

~Tz.  " - z.  1 
i+l  l-l 


i+  1 ■ i 

z.  , - z. 
i+ 1 i 


■ i i-l 

z.  - z. 
i i-l 


V-  i - v.  , 
1+1  1-1 


- A(1)  = S. 


i->- 1 


i-l 


(al 


i = 2 (N-l)  , 


where  K.--  K(z.),  v.=v(z.),  z =0,  and  z..  1, 
i i i i 1 X 

This  is  the  version  employed  in  AFGWC-BI.M.  It  is  based  on  the  notion  that 
a differentiation  at  a point  is  to  be  replaced  by  the  mean  rate  of  change  of  the 
quantity  over  the  interval  between  the  two  points,  one  on  each  side,  next  to  the  point 
of  concern.  It  is  readily  shown  that  this  formulation  is  consistent  and  that  the  local 
truncation  error  is  given  by 


fli 


= A1  - (K(z)y'). 
1 1 

ih.(v"\  + 

L * i 


+ 0(Az2)  . 


i(K  (y"l.  + A(K  (y'l. 


(zi  , - 


2 ?.  i-  /.  ,1 

l i-l 


(61 


Type  2 is  obtained  by  approximating  the  value  of  K(z)  at  the  midpoint  of  an 
interval  by  the  average  of  its  values  on  the  two  boundaries  of  the  interval.  It  is 
given  bv 


~ , - 7 7 

i+i  i-i 


K.  + K 

1 -V- 


i+ 1 


L 


Vi+  1 ~ Vi 
Zi+1  ‘ Z\ 


K.  , + K. 

i-l  i 


% i ~ vi-  1 

7,  - Vl 


, (21 


(71 


< N - 1 1 . 


It  mav  also  be  shown  that  the  Tvpe  2 formulation  is  consistent  and  has  the  same 
truncation  error  up  to  the  first  order  as  does  Tvpe  1;  that  is. 


fzi  = A(2)  - (K(z)y').  = Tu  + 0(Az2)  . (8) 

Notwithstanding  these  similarities,  there  is  a subtle  difference  between  the  two 
formulations.  The  difference  is  revealed  in  the  characteristics  of  the  resultant 
coefficient  matrices  that  become  pertinent  only  when  the  steady-state  solution  is 
sought.  It  becomes  more  obvious  when  Eqs.  (5)  and  (7)  are  rewritten,  respectively. 


z.  - z.  , ) 

l l- 1 


(Ki+l-Ki-l) 
(Zi+l  " Zi-1 


z.  , - z.)  z.  - z.  , 

1+1  1 1 1-1 


<Ki.l  -Kl-1> 


* <V.- V 4 ^1.1  - v.1  4 


(K.  + K.  ,) 

l l- 1 

(zi  ' zi-l’ 


(K.  + K.  , ) (K.  + K.  .) 

i i - 1 + ( i+l 

(zi  ' zi-l’  (Vl  ' zi' 


(K.  + K ) 

+ " (zi-i  - zi-i,si 


i = 2, (N-l) 


So  long  as  K.>  0,  for  all  i’s,  as  has  been  assumed  in  the  definition  of  the  prob- 
lem, the  coefficients  a.,  of  y,  in  the  i-th  equation  in  Eq.  d')  are  seen  to  satisfy 

aii<  °:  aik£  0 if  i^k;  and  ^ la.J  = I a. . I for  all  i=2 (N-l).  The  same,  how- 

k\  i . 

ever,  does  not  hold  true  for  the  coefficients  a.^  of  y^  in  Eq.  (5  ). 

This  means  that  when  the  system  of  equations  for  v,  (incorporating  the  boundary 
conditions)  is  obtained,  the  coefficient  matrix  formed  from  Eq.  (7,'  is  irreduciblv 
diagonally  dominant,  irrespective  of  the  K profile  and  of  the  computational  grid 
involved,  and  is,  therefore,  always  non-singular.  The  same  cannot  be  said  of  the 
corresponding  coefficient  matrix  formed  from  Eq.  (5  ).  This,  however,  does  not 
imply  that  the  latter  is  necessarily  singular.  In  fact,  with  the  restricted  structure 
that  it  has,  it  is  rather  a fortuitous  and  rare  occurrence  that  the  coefficient  matrix 
associated  with  Eq.  (5  # ) turns  out  to  be  singular  in  practice.  The  fact  that  one  is 
and  the  other  is  not  diagonally  dominant  is  obviously  sufficient  to  produce  a 
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significant  difference  in  discretization  error  between  the  two  formulations,  as 
described  below. 

Two  typical  patterns  of  discretization  error  seem  to  emerge  from  all  the  cases 
considered.  Figures  4a  and  4b  present  each  of  these  types,  using  the  results  ob- 
tained with  the  computational  grid  "10U"  that  subdivides  the  domain  1-  z-  0 into  10 
equal  subintervals.  In  each  of  these  figures  there  are  three  entries:  the  continuous 
solution  yg(z)  shown  by  a solid  curve;  the  discretization  error  with  Type  1 
finite-difference  formulation  at  each  grid  point  represented  by  O ; and  the  discre- 
tization error  with  Type  2 formulation  represented  by  x.  Note  that  while  the  range 
and  scale  of  y^  remain  the  same  in  both  figures,  those  for  discretization  error  vary 
from  one  to  the  other. 


Figure  4a.  Pattern  A of  the  Discretiza- 
tion Errors  of  the  Finite-Difference 
Formulations 


0 0 2 0.4  0.6  0.8  1.0 

0 0.02  0 04  0.06  008  01 


Y—  0 02  04 

e — - - 0 06  -004  -002 


Figure  4b.  Pattern  B of  the  Discretiza- 
tion Errors  of  the  F inite- Difference 
Formulations 


08  1 0 
0 02  0 04 


Pattern  A,  typified  by  those  shown  in  Figure  4a,  is  observed  in  eddy  diffusivi- 
ties  Kj,  and  and  in  association  with  all  the  source  functions  considered. 


The  discrete  solutions  of  both  Type  1 and  2 formulations  are  biased  toward  the 


straight  line  joining  the  two  boundary  values  yQ  and  Vj  and,  therefore,  have  pro- 
files of  great  similarity.  The  magnitudes  of  the  discretization  errors  increase 
with  the  local  curvature  of  the  continuous  solution  ys.  Differences  among  the 
eddv  diffusivities  and  among  the  source  functions  produce  differences  in  the  magni- 
tude of  individual  errors,  but  do  not  appear  to  create  am  outstanding  change 
either  in  the  distributions  or  in  the  relation  between  the  two  types. 

In  all  cases  computed,  the  error  in  Type  2 was  smaller  than  that  in  Type  1 at 
a ratio  that  remained  nearly  constant  at  all  grid  points.  Another  example,  which 
is  associated  with  a slightly  different  profile  of  continuous  solution  vw,  is  shown 
in  Figure  5a  to  illustrate  the  common  features  of  Pattern  A. 


Figure  5a.  Pattern  A of  the  Discretiza- 
tion Frrors  of  the  1 ir.ite- Difference 
1 ormulations 


Pattern  B,  characterized  bv  the  distributions  shown  in  Figure  4b,  is  observed 
in  eddv  diffusivities  Kj,  K_,  and  Kp.  One  essential  difference  of  this  type  from 
Type  A,  readily  noticeable  in  Figure  4b,  is  the  contrast  between  the  two  types  of 
formulations.  This  contrast  becomes  starker  as  eddy  diffusivity  is  changed  from 
Kj  to  Kj.,  as  showm  in  Figure  5b.  The  only  difference  betw'een  and  K,.  is  the 
minimum  value  of  eddv  diffusivitv  (min  K,  = 1,  min  K,  = 0.  1>.  As  shown  earlier 
in  Figure  1b,  there  is  little  difference  in  profile  between  and  K,.  KR,  which  has 
the  minimum  of  0.01,  is  indistinguishable  from  K,.  on  the  scale  of  Figure  lb.  How  - 
ever, when  Kp  is  used,  the  Type  1 discretization  error  becomes  so  large  that  it 
cannot  be  put  in  Figure  5b. 

(n  Tvpe  B,  the  discretization  error  of  Type  2 formulation  ri  iains  w'ell- 
behaved,  while  that  of  Tvpe  1 shows  a distinct  mark  of  pathological  disorder.  Table 
1 illustrates  this  bv  presenting  the  errors  of  the  largest  magnitude  for  both  formu- 
lations under  a variety  of  conditions.  To  be  sure,  the  discretization  errors  in  both 
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formulations  become  larger  as  the  approximate  representation  afforded  by  the 
discrete  system  becomes  worse.  It  is,  however,  clear  that  Type  1 formulation  is 
definitely  inferior  to  Type  2. 


Figure  5b.  Pattern  B of  the  Discretiza- 
tion Errors  of  the  Finite-Difference 
Formulations 


-4 

Table  1.  Errors  of  Largest  Magnitude  (X  10  ) in  Pattern  B 


Eddy  Diffusivity 

Source 

K4 

K . 

Kfi 

Function 

Type 

1 Type  2 

Type  1 

Type  2 

Type  1 

Type  2 

S1 

347 

-116 

-5863 

-721 

-44443 

1824 

S2 

273 

-89 

-8190 

560 

-60903 

2447 

S3 

-319 

-85 

-7624 

671 

-56635 

2263 

S4 

-377 

108 

-8788 

829 

-65163 

2627 

S5 

27  3 

-90 

-6760 

554 

-50535 

2041 

S6 

449 

-167 

-4102 

-950 

-32251 

1385 

A close  examination  reveals  that  the  pattern  exhibited  by  the  discretization 
error  of  Type  1 formulation  is  closely  associated  with  either  weak  or  the  lack  of 
diagonal  dominance  in  the  coefficient  matrix  of  the  system  of  equations  that  deter- 
mines the  discrete  solution.  For  example,  even  though  K4  provides  a diagonally 
dominant  coefficient  matrix,  the  degree  of  dominance  at  grid  points  near  the  level 
of  the  minimum  diffusivity  is  less  than  that  of  those  in  the  vicinity  of  the  maximum 
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diffusivity  in  Kj,  K 2>  or  Kg.  In  Kg  and  Kg  diagonal  dominance  does  not  exist  at  all 
at  the  two  grid  points  (levels  0.  3 and  0.  4)  adjacent  to  the  level  of  minimum  diffusiv- 
ity (z=0.32),  and  the  disparity  is  greater  in  Kg  than  in  Kg.  This  is  the  reason  why 
the  disorder  is  observed  in  the  discretization  error  of  the  Type  1 formulation  and 
not  in  that  of  Type  2,  which  preserves  diagonal  dominance  under  all  conditions  by 
nature  of  its  construction.  We  may  also  note  in  Table  1 that  the  error  of  the  largest 
magnitude  for  the  Type  2 formulation  is  inevery  case  smaller  than  that  for  Type  1. 

3.2  Structure  of  Computational  Grids 

Given  the  domain  and  range  of  prediction,  density  of  the  computational  grid  for 
a dynamical  prediction  model  is  restricted  within  a well-defined  range  that  is  largely 
determined  by  such  factors  as  size  and  speed  of  the  computer  employed  and  the 
preparation  time  available  for  the  required  forecast.  With  the  understanding  that 
AFGWC-BLM  was  designed  to  make  maximum  use  of  the  computer  capability,  we  com- 
pared three  computational  grids  for  accuracy  using  the  Type  2 finite -difference  for- 
mulation. They  are:  "10U",  the  grid  that  subdivides  the  domain,  l-z^0,  into  10 
equal  subintervals;  "7U",  the  grid  that  subdivides  the  domain  into  7 equal  subinter- 
vals; and  "7N",  the  grid  consisting  of  levels  z = 1/30,  1/10,  1/5,  2/5,  3/5,  4/5, 
and  1.  The  last  grid  very  closely  resembles  the  vertical  grid  of  AFGWC-BLM,  and 
is  intended  to  represent  the  latter  in  relation  with  the  other  two. 

Figure  6 assembles  six  different  samples  selected  to  represent  the  whole  spec- 
trum of  variation  within  the  total  population  considered.  The  discretization  error 
at  each  grid  point  of  the  steady-state  solution  obtained  with  each  of  the  three  grids 
is  shown  in  these  samples.  It  is  obvious  that  the  relative  merit  of  each  of  the  three 
grids  changes  from  one  case  to  the  next.  It  is,  however,  not  too  difficult,  even  on 
the  basis  of  these  six  samples,  to  arrive  at  an  overall  judgment  that  ranks  "10U", 
"7U",  and  "7N"  in  the  order  of  decreasing  accuracy. 

Table  2 summarizes  the  mean  absolute  errors  of  the  three  grids  for  six  diffu- 
sivities  and  three  source  functions.  Since  the  errors  are  defined  on  a different  set 
of  grid  points  for  a different  grid,  we  calculated  various  statistics  and  tried  to  see 
if  a meaningful  comparison  could  be  made  with  such  statistics.  Thus,  in  each  eddy 
diffusivity  in  Table  2,  the  first  line  gives  the  nine -point  average  of  absolute  errors 
from  "10U",  the  second  the  six -point  average  from  "7N",  and  the  third  the  six -point 
average  from  "7U".  Because  of  the  differences  in  composition  of  these  averages, 
they  are  not  directly  comparable.  However,  the  three  values  on  the  last  line  (the 
first  two  being  the  mean  absolute  errors  from  z = 1/5,  2/5,  3/5,  and  4/5  that  are 
common  to  both  "10U"  and  "7N",  and  the  last  the  average  of  absolute  errors  from 
z = 2/7,  3/7,  4/7,  and  5/7  on  "7U")  appear  to  be  quite  compatible  with  each  other, 
particularly  when  taken  with  cognizance  of  the  other  line  in  each  column.  The  table 
strongly  suggests  that  "10U”  can  reduce  the  discretization  error  measured  by  mean 
absolute  error  to  half  that  of  "7N",  and  even  "7U"  can  be  reduced  bv  approximately 
20  percent  in  the  majority  of  cases  considered. 
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Source  Function 


Eddy 

Diffusivity 

S1 

10U  7N  7U 

S4 

10U  7N  7 U 

S6 

10U  7 N 7 U 

K1 

29829 

37378 

41143 

29121  52727  39534 

39430 

49730 

54668 

38564  70144  52635 

26349 

42866 

36652 

25387  61479  35347 

K2 

25977 

33873 

36459 

25660  47006  35655 

3827  1 

49913 

54151 

37856  69617  53127 

21574 

27568 

30637 

21351  38691  30128 

K3 

11215 

21676 

15413 

9694  30488  9615 

45807 

63898 

65193 

46174  82227  65795 

20118 

97  50 

26861 

19661  10733  26785 

K4 

5203 

18503 

10816 

6933  22084  13011 

4649 

15603 

9138 

4948  17859  9643 

5450 

20823 

14586 

8302  25799  18933 

K5 

26640 

53249 

43915 

34395  69404  55702 

22693 

44927 

46661 

28666  57104  59137 

27  660 

53399 

45200 

41181  80834  57495  ' 

K6 

56457 

73500 

77  193 

57454  97266  101899 

60660 

63794 

57261 

49729  82966  7 1249 

57835 

83804 

79874 

67705  112428  105100 

So  far,  our  attention  has  been  on  the  synopsis  of  discretization  errors  over  the 
entire  domain.  Since  the  nonuniformity  of  "7N"  arises  from  having  z = 1/30  and 
z * 1/10  in  the  computational  grid,  one  might  ask  what  such  inclusion  would  do  to 
the  accuracy  of  the  local  estimates  at  these  levels.  Comparisons  of  discretization 
errors  were,  therefore,  made  among  the  three  computational  grids  by  estimating 
y at  z = 1/30  and  1/10  from  the  values  at  z = 0 and  1/7  through  the  linear  inter- 
polation on  "7U",  and  by  estimating  y at  z - 1/30  from  the  values  at  z = 0 and  1/10 
through  the  linear  interpolation  on  "10U". 

Table  3,  in  which  the  values  of  the  continuous  solution  and  the  errors  of  esti- 
mate in  three  grids  are  given,  summarizes  the  results  for  six  diffusivities  and 
three  source  functions.  It  is  found  that  the  estimates  in  "10U"  are  more  accurate 
than  those  in  others  at  both  levels  in  the  great  majority  of  the  cases.  It  is  also  true 
that  the  estimate  in  "7U"  at  z - 1/10  is  generally  better  than  that  in  "7N".  On  the 
other  hand,  at  z = 1/30,  "7N"  and  "7U"  are  about  even  in  accuracy. 


18 


- fi 

Table  3.  Local  Discretization  Errors  at  z = 1/30  and  1/10  (X  10  ) 


I 

Putting  these  results  together,  it  is  clear  that  we  have  not  only  failed  to  uncover 
any  definite  advantage  of  the  nonuniform  grid,  but  also  have  convincingly  shown  that 
a uniform  grid  is  preferred  to  a nonuniform  one,  even  at  a density  as  low  as  those 
considered  in  the  present  investigation,  unless  there  is  a definite  and  stringent  re- 
striction on  the  parameter  functions  that  makes  the  choice  of  an  uneven  grid  in- 
disputably more  preferrable. 
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Source  Function 

S 

s 

s„ 

Eddy 

1 

4 

6 

Diffusivity 

Level 

10U 

7N 

7 U 

10U 

7 N 

7 U 

10U 

7 N 

7 U 

y 

17839 

6287 

24089 

K1 

1/30 

e 

1320 

3394 

1954 

3138 

4516 

4486 

1034 

2849 

1334 

1/10 

y 

52183 

20655 

70259 

e 

5694 

9968 

7 195 

7619 

13288 

11663 

5110 

8430 

6009 

y 

26539 

6203 

37475 

1/30 

K„ 

e 

132 

487  2 

-79 

3893 

5503 

5499 

-1005 

2771 

-1836 

2 

1/10 

y 

74097 

21236 

104274 

e 

5917 

10341 

5283 

9052 

15507 

13871 

4986 

787  6 

2642 

y 

128936 

17  97 

195136 

1/30 

K- 

e 

-55876 

1589 

-68538 

11981 

18743 

14384 

-87755 

-7  287 

-107527 

3 

1/10 

y 

240957 

18820 

363361 

e 

-21777 

6514 

-59762 

22513 

35737 

29722 

-41217 

-8282 

-10053 

y 

12164 

-3045 

18536 

1/30 

K, 

e 

4346 

4965 

8331 

3748 

4738 

7497 

5598 

47  94 

10442 

4 

1/10 

y 

45955 

-2865 

69207 

e 

3574 

177  16 

15529 

4973 

17444 

16220 

3194 

16945 

17726 

y 

3793 

-1237  1 

9961 

1/30 

Ks 

e 

5349 

87  60 

10041 

4236 

8451 

8330 

707  5 

8835 

12973 

5 

1/10 

y 

14770 

-38241 

37840 

e 

12656 

33115 

26732 

13836 

32698 

26118 

13268 

33212 

30962 

y 

1173 

-15030 

7 123 

1/30 

K 

e 

6801 

10782 

11688 

5527 

10406 

97  81 

8774 

11063 

14913 

6 

1/10 

y 

4582 

-48647 

26873 

e 

19340 

41155 

34002 

20139 

40496 

32899 

20819 

42046 

39234 
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4.  CONCLUSIONS 

Using  discretization  error  in  the  steady -state  solution  as  a measure  of  in- 
accuracy for  the  time -dependent  solution,  we  investigated  the  relative  accuracy  of 
discrete  solutions  obtained  with  a number  of  discretization  procedures  that  were 
comparable  in  operational  aspects  with  that  employed  in  AFGWC-BLM.  We  have 
come  to  the  following  two  major  conclusions: 

(1)  The  finite -difference  formulation  presently  used  in  AFGWC-BLM  has  an 
undesirable  feature  that  can  be  removed  entirely  by  another  formulation.  The 
latter  is  of  the  same  accuracy  in  local  truncation  error  and  requires  about  the 
same  computational  effort  as  the  present  one. 

(2)  The  computational  grid  currently  employed  in  AFGWC-BLM  not  only  shows 
no  distinct  advantage,  it  also  exhibits  characteristics  that  are  considered  to  be 

S inferior  to  those  of  the  uniformly  space  grid.  Both  types  of  grids  require  the  same 

I amount  of  storage  space  and  about  the  same  amount  of  computation  time. 

’ 


